On the origin of the Fermi arc phenomena in the underdoped cuprates: signature of 

KT-type superconducting transition 
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We study the effect of thermal phase fluctuation on the electron spectral function A{k, uj) in a 
d-wave superconductor with Monte Carlo simulation. The phase degree of freedom is modeled by 
a XY-type model with build-in d-wave character. We find a ridge-like structure emerges abruptly 
on the underlying Fermi surface in A{k, uj = 0) above the KT-transition temperature of the XY 
model. Such a ridge-like structure, which shares the same characters with the Fermi arc observed 
in the pseudogap phase of the underdoped cuprates, is found to be caused by the vortex-like phase 
fluctuation of the XY model. 
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The origin of the pseudogap phenomena in the under- 
doped cuprates is among the most hotly debated issues 
in the high-Tc physics. Existing theories on the pseudo- 
gap falls into two categories: while it is taken by many 
as the evidence for the existence of preformed Cooper 
airs and strong phase fluctuation in the normal state jl!- 
_^ some people believed that it should be attributed to 
a yet another unidentified order j^Q. One puzzle about 
the pseudogap is its momentum and temperature depen- 
dence. While early experiments suggests that the pseu- 
dogap may inherit the d-wave structure of the supercon- 
ducting gap, more detailed ARPES measurement show 
that this is not the case^-Q- Below T* , the pseudo- 
gap is first seen in a small momentum region around 
the antinodal point (tt, 0). The gaped region enlarges 
with decreasing temperature until the whole Fermi sur- 
face(except the nodal points) is gaped below Tc- The 
gradual development of the pseudogap on the Fermi sur- 
face leaves the system with a finite segment of ungapped 
Fermi surface which is called Fermi arc in the literature. 
The length of the Fermi arc increases with temperature 
and is reported to have a jump at T^: the Fermi arc 
emerges abruptly at T^. 

A metal with a unclosed Fermi surface is a strange 
animal in the zoo of the Landau Fermi liquid theory. The 
Fermi surface, by its very definition as an equal-energy 
contour at the chemical potential, is always closed for any 
well defined quasiparticle dispersion. It is proposed that 
the Fermi arc is just a half of a closed pocket-like small 
Fermi surface whose other half is too weak to be detected 
by ARPES [10]. Such a scenario is widely adopted in the 
second kind of theories for the pseudogap (related to an 
unidentified order). In this paper, we try to understand 
the Fermi arc phenomena in the framework of the phase 
fluctuation scenario. In this scenario, the Fermi arc is not 
a genuine Fermi surface in the momentum space across 
which a finite jump of the occupation number occur, but 
a phase fluctuation induced pile up of the spectral weight 
on the Fermi energy in a d-wave paired state. 

The effect of the phase fluctuation on the electronic 



spectral function have been studied by many authors [lll- 
lj|. However, most of these studies are limited to the 
semiclassical treatment of the phase fluctuation. At the 
semiclassical level, the phase fluctuation is modeled by a 
uniform supercurrent with a phenomenological distribu- 
tion of velocity. The Doppler effect of the quasiparticle 
moving in such a uniform supercurrent background then 
causes electronic spectral function to shift in both the 
momentum and energy. More specifically, when the cen- 
ter of mass momentum of the moving Cooper pairs is 
2q, the electron pairing then happens between momen- 
tum q + k and q — k, and the excitation energy become 
Ek — Vk ■ q to first order in q. The net effect of the su- 
percurrent on the electron spectral function is to shift it 
by q in momentum and w/c • g in energy. Here Vk is veloc- 
ity of the electron derived from the bare band dispersion 
relation. 

In a recent study [l3|, the author made the important 
observation that in a d-wave superconductor the Doppler 
effect of the fluctuating supercurrent will cause pile up 
of the spectral weight exactly at the Fermi energy on the 
underlying Fermi surface, no matter the detailed distri- 
bution of the fluctuating supercurrent. The pile up of the 
spectral weight results in inverse square root singularity 
at the Fermi energy in the electronic spectral function, 
resembling a quasiparticle peak. The key point for this 
to happen is that in a d-wave superconductor, the quasi- 
particle has linear dispersion around the gap node. For 
such a system, the effect of the momentum shift and the 
energy shift caused by the fluctuating supercurrent can- 
cels each other for q in the nodal direction. For q in 
general direction, the transverse component of q±^ causes 
no energy shift in leading order and simply move the gap 
node along the underlying Fermi surface transversely by 
a distance of q±. 

The semiclassical approach, though intuitively attrac- 
tive, has limited value to give quantitative estimate of 
the phase fluctuation effect. To simulate the actual tem- 
perature dependence of the phase fluctuation effect, one 
must go beyond the long wave length limit and take into 
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account of the singular effect of the vortex excitation on 
the quasiparticle motion. Even for the simulation of the 
phase fluctuation involved in the far field of the vortex, 
the semiclassical approach has the problem that it is un- 
able to provide a quantitative measure of the strength of 
the fluctuation. 

In this paper, we report the results of the Monte Carlo 
simulation of the effect of thermal phase fluctuation on 
the electronic spectral function in d-wave superconduc- 
tors. The phase fluctuation is described in our work by 
a XY-type model with a build-in d-wave character. The 
phase fluctuation then couples to the quasiparticles of the 
d-wave superconductor via a standard pairing model with 
the magnitude of the pairing potential A fixed. The sim- 
ulation shows clearly how the Fermi arc emerge abruptly 
above the KT transition temperature of the XY model 
and then increase in strength with temperature. These 
results indicate that the Fermi arc observed in the under- 
doped cuprates is a phase fluctuation effect. The abrupt 
emergence of the Fermi arc above Tc suggests that the 
superconducting transition of the underdoped cuprates 
are of KT transition-type. 

The phase degree of freedom is described by the fol- 
lowing XY model 



H. 



phase 



= J ^ cos( 
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in which 9a is the phase variable defined on the bonds of 
a square lattice on which the electrons reside. The center 
of these bonds form another square lattice(see Fig.l) and 
the sum in Eq.(l) is over neighboring sites of the latter 
square lattice. The coupling constant J is set to be posi- 
tive so that a d- wave-like phase structure is favored at low 
temperature in this model, namely the phase variable on 
the neighboring bonds tends to have a phase difference 
of TT. Similar models are also used in other works 12, 15| . 

In principle, the energetics of the phase degree of free- 
dom should be determined by the dynamics of the elec- 
tron system itself. Here we have treated them as sep- 
arated degrees of freedom governed by a coupling con- 
stant J. Except for being more simple to simulate, such 
a model has the advantage that it takes into account the 
fact that in the underdoped cuprates the spin pairing 
and phase coherence are controlled by two separate en- 
ergy scaled, [^. Such a separation of the energy scales is 
believed to the key for the formation of the pseudogap 
phase in most scenarios. 

The electron degree of freedom is described by the 
standard pairing model of the form Hp^ir = Ht -|- Ha, 
in which Hf and Ha are given by 



<i,j>,(J 




FIG. 1: The lattice on which model Eq.(l) is defined on. The 
electrons reside on the square lattice made of the solid lines. 
The phase variables resides on the bonds of this square lattice 
which are indicated by the gray dots. These gray dots form 
another square lattice of itself. 
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in which t and t' are the hopping integral of the elec- 
tron on the square lattice between neighboring and next 
neighboring sites. The next neighboring hopping t' is 
introduced to model the curvature of the real Fermi sur- 
face. A denotes the magnitude of the pairing potential 
and is kept as a constant in the whole simulation. 9 a is 
the phase angle on the bond between site i and j. 

In the simulation, we first generate a thermal phase 
fluctuation configuration from the distribution p — 
e:>q3{—'51phase/kBT). Then we diagonalize the pairing 
Hamiltonian Eq.(2) with the generated phase fluctuation 
configuration and calculated electronic spectral function. 
We then average the electronic spectral function on dif- 
ferent phase fluctuation configurations. 

For each individual phase fluctuation configuration, 
the translational invariance of the pairing Hamiltonian is 
lost. However, when the result is averaged over all pos- 
sible phase fluctuation configurations the translational 
invariance is recovered. In the simulation, we have av- 
eraged the electron spectral function over all configura- 
tions generated from a given one by space translation at 
each step. After the average translational invariance is 
recovered and electronic spectral function as a function 
of momentum A{k, io) is given by 



A(k,uj) 



(3) 



in which 



denotes the eigenvectors 



of the following B-dG equation with eigenvalue -E", 



«l,J>>,cr 
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in which Ht^i^ = -tStj+R^ - t'Sij+R^, - fiSij and 
HA,i.j — Ae'^°5ij_|_/j^ are the matrix element of Ht and 
Ha- Here sum over repeated indexes is assumed and Rs 
and Rs' denotes the neighboring and next neighboring 
vector on the square lattice. When A{k,uj) is calculated 
for a given phase fluctuation configuration, we do the 
average over the configurations by the standard Monte 
Carlo procedure 



{Aik,u)) 
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A{k,Lu) 



(5) 
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FIG. 2: The correlation(defined in the text) between the 
phase variables on the most far apart bonds of a 64 x 64 
lattice as a function of temperature. The abrupt drop of this 
correlation around ksT/J ~ f indicates the KT transition of 
that finite lattice. 

The diagonalization of the B-dG Hamiltonian is the 
most time consuming part of our simulation. At the same 
time, to resolve the Fermi surface on the finite momen- 
tum mesh, we have to study lattice of large enough size. 
To study relatively larger lattice so that the Fermi surface 
can be resolved, we adopt the following cutoff strategy. 
First, we Fourier transform(by FFT algorithm) the B- 
dG Hamiltonian into momentum space. The transformed 
Hamiltonian reads 



-^k5k,k' 



(6) 



At' cos kx cos ky — fi, 



in which = — 2t(cos kx + cos ky 
Ak,k' = l/NJ2<^,J>^e'''^'e''''^^-'''''^'. To simulate the 
Fermi arc phenomena, we only need to calculate the 
electron spectral function right on the Fermi energy, 
A{k,io = 0). Obviously, those momentums for which 
^ A are not expected to contribute significantly to 
A{k,uj — 0). For this reason, we set a cutoff energy Ec 
and neglect all momentums for which > Ec- We then 
diagonalize the truncated B-dG Hamiltonian below the 
cutoff energy and calculate the electron spectral function 
on the Fermi energy. The convergence of such a cutoff 



strategy can be easily checked by varying the cutoff en- 
ergy Ec- In real simulation, we find Ec = 5 A is large 
enough. 

We now present the results of our simulation. The 
simulation is done on a 64 x 64 lattice. We have used 
t'/t = —0.3 and A/t = 0.1 to model the real system. 
The coupling strength of the phase variables, J, is set 
as the unit for temperature. The strength of the phase 
fluctuation can be read out from the correlation function 
of the phase variables. In Figure 2, we plot the corre- 
lation C{T) = {cos{9a — Oa')) between the phase vari- 
ables on the most far apart bonds of the 64 x 64 lattice 
as a function of temperature. For infinite system, such 
correlation is zero at finite temperature as required by 
the Wagner-Mermin theorem. The abrupt drop of this 
correlation signifies the KT transition of the finite sys- 
tem. For our 64 x 64 lattice, we find the abrupt drop of 
the phase correlation occurs around ksT/ J ^ I, slightly 
above the KT transition temperature for the infinite sys- 
tem {ksT ~ 0.8923J)[il,[i3]. 

To achieve statistical independence of the samples gen- 
erated from the Monte Carlo procedure, we have used 
10^ local update of the phase variables to generate every 
new phase fluctuation conflguration. The flrst 10* con- 
figurations are discarded for thermalization. The elec- 
tron spectral function is averaged over 1000 indepen- 
dent phase fluctuation configurations. Figure 3 shows 
the A{k,uj = 0) calculated at six temperatures. Among 
the six temperatures, one{kBT = 0.5J) is below the KT 
transition point and the other five are above itlksT ~ 
J, I.IJ, 1.2J,1.3J, 1.4J). 

Below the KT transition temperature(rR-T), the phase 
fluctuation effect is almost totally quenched and the elec- 
tron spectral function at the Fermi energy is composed 
of four sharp peaks at the four nodal points, as in an 
ideal d-wave superconductor. Right above the KT tran- 
sition temperature, a ridge-like structure in A(k,uj = 0) 
emerges on the underlying Fermi surface The inten- 
sity of this ridge-like structure increases rapidly in a small 
temperature region between fc^T = J and ksT = I.IJ. 
For ksT > I.IJ, the increase of electron spectral func- 
tion on the Fermi energy becomes much slower. Such 
temperature dependence of the electron spectral func- 
tion at the Fermi energy is in close resemblance to that 
observed in experiments @. It should be noted that the 
ridge-like structure emerges simultaneously on the whole 
Fermi surface, though it is clear anisotropic on the Fermi 
surface. However, considering the finite resolution of the 
ARPES measurement, it is still reasonable to deflne a 
finite segment of the Fermi surface as the Fermi arc for 
each temperature by setting some threshold for the spec- 
tral weight. The arc length defined in such a way clearly 
increases with temperature. 

To quantify the temperature dependence of the elec- 
tron spectral function in a more objective way, we plot in 
Figure 4 the integration of the electron spectral function 
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FIG. 3: The temperature dependence of the electron spectral function at the Fermi energy A{k, cj — 0). (a)fcsr = 0.5 J, 
(b)A;sT = J, {c)kBT = l.lj, (d)fcflr = 1.2J, (e)itflT = 1.3J,(f)fei3r = 1.4J. Shown here are the results in the first quadrant 
of the Brillouin zone only. 

at Fermi energy in the Brillouin zone W{io = 0,r) = 
J dkA{k, = 0) as a function of temperature. The 
quench of the phase fluctuation effect below Tkt and 
the abrupt emergence of the ridge-like structure above it 
can be clearly seen in this figure. 

We thus conclude that the main features of the exper- 
imental observations concerning the Fermi arc phenom- 
ena can be accounted for by our phase fluctuation model. 
The quench of the phase fluctuation effect below Tkt and 
the its abrupt emergence above Tkt indicates that Fermi 
arc phenomena is controlled by the same mechanism that 
induce the KT transition, namely the vortex excitation. 
Indeed, according to the study of the XY model, the 
vortex density enjoys a proliferation just above the KT 
transition point. With these understanding, it is tempt- 
ing to argue that the Fermi arc behavior observed in the 
high-Tc cuprates provide strong evidence that the super- 
conducting transition in the underdoped cuprates are of 
KT transition-type. 
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FIG. 4: The integrated electron spectral function over mo- 
mentum at the Fermi energy W{uj = 0,T) = J dkA{k, lj = 0) 
as a function of temperature. 
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